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(Dated: August 16, 2006) 

We present the first experimental discrimination between the 2e/3 and 4e/3 top quark electric 
charge scenarios, using top quark pairs (tt) produced in pp collisions at y / s=1.96 TeV by the Fermilab 
Tevatron collider. We use 370 pb _1 of data collected by the DO experiment and select events with 
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at least one high transverse momentum electron or muon, high transverse energy imbalance, and 
four or more jets. We discriminate between 6- and 6-quark jets by using the charge and momenta 
of tracks within the jet cones. The data is consistent with the expected electric charge, |g| — 2e/3. 
We exclude, at the 92% C.L., that the sample is solely due to the production of exotic quark pairs 
QQ with \q\ — 4e/3. We place an upper limit on the fraction of QQ pairs p < 0.80 at the 90% C.L. 

PACS numbers: 13.85.Rm, 14.65.Ha 



The heavy particle discovered by the CDF and DO col- 
laborations at the Fermilab Tevatron proton-antiproton 
collider in 1995 [1] is widely recognized to be the top 
quark. Currently measured properties of the particle are 
consistent with standard model (SM) expectations for 
the top quark. However, many of the properties of the 
particle are still poorly known. In particular, its electric 
charge, a fundamental quantity characterizing a particle, 
has not yet been determined. 

To date, it is possible to interpret the discovered par- 
ticle as either a charge 2e/3 or — 4e/3 quark. In the 
published top quark analyses of the CDF and DO col- 
laborations [2|, there is a two- fold ambiguity in pair- 
ing the 6-quarks and the W bosons in the reaction 
pp — > ti — > W + W~bb, and equivalently, in the electric 
charge assignment of the measured particle. In addition 
to the SM assignment, t — > W + b, "i" — » W~b is also con- 
ceivable, in which case "f" would actually be an exotic 
quark, Q, with charge q = — 4e/3 (charge-conjugate pro- 
cesses are implied). It is possible to fit Z — > £ + £~ and 
Z — > 66 data assuming a top quark mass of m t — 270 GeV 
and a right-handed 6-quark that mixes with the isospin 
+1/2 component of an exotic doublet of charge — le/3 
and — 4e/3 quarks, (Ql ,Q4)r 0]. In this scenario, the 
— 4e/3 charge quark is the particle discovered at the Teva- 
tron, and the top quark, with mass of 270 GeV, would 
have so far escaped detection. 

In this Letter, we report the first experimental discrim- 
ination between the 2e/3 and 4e/3 charge scenarios. We 
also consider the case where the analyzed sample con- 
tains an admixture of SM top quarks and exotic quarks 
and place an upper limit on the exotic quark fraction. 
Our search strategy assumes each quark decays 100% of 
the time to a W boson and a 6-quark. We use the lepton- 
plus-jets channel which arises when one W boson decays 
leptonically and one decays hadronically. The charged 
leptons (e//i) originate from a direct W decay or from 
W —* r — > e/fi. We require that the final state have at 
least two 6-quark jets. The data used in this Letter were 
collected by the DO experiment from June 2002 through 
August 2004 and correspond to an integrated luminosity 
of 370 pb" 1 . 

The DO detector includes a tracking system, calorime- 
ters, and a muon spectrometer The tracking sys- 
tem is made up of a silicon microstrip tracker (SMT) 
and a central fiber tracker, located inside a 2 T super- 
conducting solenoid. The SMT, with a typical strip 
pitch of 50-80 /xm, allows a precise determination of 



the primary interaction vertex (PV) and an accurate 
determination of the impact parameter of a track rela- 
tive to the PV [5J]. The tracker design provides efficient 
charged-particle measurements in the pseudorapidity re- 
gion 1 77 1 < 3 [(|. The calorimeter consists of a barrel 
section covering |^| < 1.1, and two end caps extending 
to 1 77 w 4.2. The muon spectrometer encapsulates the 
calorimeter up to |?7| = 2.0 and consists of three layers of 
drift chambers and two or three layers of scintillators 0] • 
A 1.8 T iron toroidal magnet is located outside the in- 
nermost layer of the muon detector. 

We select data samples in the electron and muon chan- 
nels by requiring an electron with transverse momentum 
Pt > 20 GeV and < 1.1, or a muon withp-r > 20 GeV 
and I77I < 2.0. The leptons are required to be isolated 
from other particles using calorimeter and tracking in- 
formation. More details on the lepton identification and 
trigger requirements are given in Ref. [8j . W boson candi- 
date events are then selected in both channels by requir- 
ing missing transverse energy, $t, in excess of 20 GeV 
due to the neutrino. To remove multijet background, 
fix is required to be non-collinear with the lepton di- 
rection in the transverse plane. Jets are defined using 
a cone algorithm [j| with radius A1Z = 0.5 [l(J. These 
events must be accompanied by four or more jets with 
Pt > 15 GeV and rapidity \y\ < 2.5. After all the above 
selection requirements are applied, we have a total of 231 
(277) events in the muon (electron) channel. 

We use a secondary vertex tagging (SVT) algorithm 
to reconstruct displaced vertices produced by the decay 
of B hadrons. Secondary vertices are reconstructed from 
two or more tracks satisfying: px > 1 GeV, > 1 hits 
in the SMT layers, and impact parameter significance 
d C a/&d ca > 3.5. A jet is considered as SVT-tagged if 
it contains a secondary vertex with a decay length sig- 
nificance Lxyjcri, > 7 11]. The determination of the 



sample composition relies on 6-tagging, c-tagging, and 
light flavour tagging efficiencies and uses the method de- 
scribed in Ref. [12]. In order to increase the purity of 
the sample we select only events with two or more SVT- 
tagged jets. In the selected sample of 21 events with two 
SVT-tagged jets, the largest (second largest) background 
is W66 (single top quark [ijj]) production with a contri- 
bution of « 5% (« 1%) to the number of selected events. 

The top or anti-top quark whose W boson decays 
leptonically (hadronically) is refered to as the leptonic 
(hadronic) top and the associated 6-quark is denoted bi 
(6/j). To compute the top quark charge we need to i) de- 
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cide which of the two SVT-tagged jets are be and bh and 
ii) determine if be and bh are b- or 6-quarks. The detected 
final state partons in the tt candidate events comprise the 
be and bh quarks, two quarks from the hadronically de- 
caying W boson, and one muon or one electron. The four 
highest-pT jets can be assigned to the set of final state 
quarks according to many permutations and there are at 
least two ways to assign the SVT-tagged jets to be and 
bh- For each permutation, the measured four- vectors of 
the jets and lepton are fitted to the it event hypothe- 
sis, taking into account the experimental resolutions and 
constraining the mass of two W bosons to its measured 
value and the top quark mass to 175 GeV. We decide 
which of the SVT-tagged jets are be and bh by selecting 
the permutation with the highest probability of arising 
from a ti event. Studies on simulated tt show that this 
gives the correct assignment in about 84% of the events. 

We measure the absolute value of the top quark charge 
on each side of the event, given by Q\ — \qe + % t | on the 
leptonic side and Q2 = \ — qe + Qb h | on the hadronic side. 
The charge of the lepton is indicated by qe, and q^, and 
qb h are the charges of the SVT-tagged jets on the lep- 
tonic and hadronic side of the event. The charges qb e 
and qi, h are determined by combining the pr and charge 
of the tracks contained within a cone of A1Z=0.5 around 
the SVT-tagged jet axis. Based on an optimization us- 
ing simulated tt events generated with ALPGEN 14\ and 
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151 ] for a full DO detector simulation, we define an 
estimator for jet charge q jct = ft ' Prf) I (J2iPr^) 
where the subscript i runs over all tracks with px > 
0.5 GeV and within 0.1 cm of the PV in the direction 
parallel to the beam axis. 

To determine the expected distributions for the top 
quark charges Q\ and Q2, it is crucial to determine the 
expected distributions for gj et in the case of a 6-quark 
or a 6-quark jet. In «5% of the ti events, one of the 
SVT-tagged jets is actually a c-quark jet arising from 
W — > cs (or its charge conjugate). Therefore we also 
need to determine the expected distribution for q^ ct in 
the case of c- and c-quark jets. 

We derive the expected distributions of jet charge from 
dijet collider data, enhanced in heavy flavor (b and c). 
We select events with exactly two jets, both SVT-tagged, 
with px > 15 GeV and \y\ <2.5. The method requires 
that the two jets are of charge conjugate flavors. To 
ensure this, we enhance bb and cc produced by flavor cre- 
ation E3, EH, by requiring the azimuthal distance 
between the jets to be larger than 3.0 and one jet (des- 
ignated as ji ) to contain a muon with px > 4 GeV. We 
refer to this sample as the "tight dijet sample," to j\ as 
the "tag jet" and to the second jet ji as the "probe jet." 

The fraction of cc events in the tight dijet sample is 
estimated using the distribution of the muon transverse 
momentum with respect to the tag jet axis (p^ 1 )- We 
fit the p r rp l distribution with a sum of two p™ 1 templates, 
one for 6-quark jets (including both prompt and cascade 



decays) and one for semi-muonic decays inside c-quark 
jets. This leads to a fraction x c of cc events of 
in the tight dijet sample and since the light flavor tag- 
ging efficiency is «15 times lower, we also conclude that 
the fraction of lighter flavor jets is negligible. The muon 
inside the tag jet comes either i) from a direct B meson 
decay, ii) a B — > D meson cascade decay, iii) an oscillated 
neutral B meson, or iv) a direct D meson decay. We find 
that further contribution from indirect D meson decay 
can be neglected. Charge flipping processes ii) and iii) 
lead to a muon of opposite charge to that of the quark 
initiating the tag jet and therefore of same sign as the 
quark initiating the probe jet. We find, with PYTHIA [n| 
simulated events and evtgen [2(| for heavy flavor de- 
cays, that charge flipping processes are x = (30 ± 1)% 
of the bb events in the tight dijet sample. This fraction 
is experimentally confirmed by studying charge correla- 
tion between muons in back-to-back muon-tagged dijet 
events. 

We denote the charge distributions for the probe jet 
when the muon on the tag side is positive or negative as 
P^+ and P„- . Similarly we define Pf to be the charge 
distribution when the jet is of flavor / = 6, 6, c, c. Given 
the fractions of cc events and of charge flipping processes 
we can write 



0.69P b 
0.30P 6 



0.30Pj, 
0.69P 5 



0.01P e 
0.01P C 



(1) 



P M + and P„- are distributions observed in data and are 
admixtures of the quark charge distributions. Equa- 
tions [T] are not sufficient to extract the four probabil- 
ity density functions (p.d.f.'s) Pf. Therefore we define 
a "loose dijet sample," where j\ is not required to be 
SVT-tagged. Using the same techniques as for the tight 
dijet sample, we find that x c — (19 ± 2)% and the same 
fraction of charge flipping processes as for the tight dijet 
sample. We refer to P^ + (P^-) as the observed p.d.f.'s 
for gj e t on the probe jet in the loose dijet sample, when 
the tag muon is positive (negative). Thus we can write 



P' 

p' 



0.567P, 
0.243P, 



0.243P 5 
0.567P5 



0.19P 5 
0.19P C . 



(2) 



We solve Eqs. 1 and 2 to obtain the P/ for 6-, 6-, c-, and 
c-quark jets. 

The P/'s are dependent on the jet pt, since pr cor- 
relates with track multiplicity in the jet, and on the 
jet y, since the tracking efficiency is rapidity-dependent. 
Therefore we must account for the different jet px and y 
spectra between the probe jets of the dijet samples and 
the 6-quark jets in preselected ti events. The P/'s ob- 
tained above are corrected by weighting the data events 
to the pt and y spectra of SVT-tagged jets in tt events. 
Figure 1(a) shows the resulting P& and P^. 

We derive the expected distributions for Qi and Q2 
by applying the assignment procedure between the SVT- 
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tagged jets and the bh, be quarks on simulated it events 
using our calculated P/'s. The true flavor / of the SVT- 
tagged jets is determined from the simulation informa- 
tion. The values of q^ h and are obtained by randomly 
sampling the distribution of Pf for the corresponding fla- 
vors. About 1% of tt candidate events contain a SVT- 
tagged light-flavor jet. In this case the p.d.f. for qj et 
is taken from simulation. In the case of a \q\ — 4e/3 
exotic quark, the expected distributions of exotic quark 
charge are derived by computing Q\ = \ — qi + qb e \ and 
Q2 = \qe + Qb h \, following the same procedure as for the 
SM top quark. The uncertainty on the mass of the top 
quark [21] is propagated as a systematic uncertainty. 

The expected distributions of Q\ and Q2 for the back- 
ground are obtained by i) performing the assignment pro- 
cedure between SVT-tagged jets and the bh, be quarks on 
Wbb simulated events, ii) using the true jet flavors / to 
sample the corresponding -P/'s. The resulting distribu- 
tions of Qi and Q2 for the background are added to the 
top charge distributions in the SM and exotic cases. We 
denote Psm (Pex) the p.d.f. 's for Qi and Q2 including the 
background contributions in the SM (exotic) case. 

For 16 of the 21 selected lepton-plus-jet events, the 
kinematic fit converges and we can assign the SVT- 
tagged jets to the be and bh quarks, thus providing 32 
measurements of the top quark charge. Figurefljb) shows 
the 32 observed values of Q\ and Q2 overlaid with the 
SM and exotic charge distributions. 

To discriminate between the SM and the exotic hy- 
potheses, we form the ratio of the likelihood of the ob- 
served set of charges qi arising from a SM top quark to 
the likelihood for the set of qi arising from the exotic 
scenario, A = [fT* PsM(ft)] / Hi p ex(<li)]- The subscript 
i runs over all 32 available measurements. The value of 
the ratio is determined in data and compared with the 
expected distributions for A in the SM and exotic scenar- 
ios. We find that the observed set of charges agrees well 
with those of a SM top quark. The probability of our ob- 
servation is 7.8% in the case where the selected sample 
contains only exotic quarks with charge \q\ — 4e/3, in- 
cluding systematic uncertainties. Thus, we exclude at the 
92.2% C.L. that the selected data set is solely composed 
of an exotic quark with \q\ = 4e/3. The corresponding 
expected C.L. is 91.2%. Table U summarizes the domi- 
nant systematic uncertainties and their cumulative effect 
on the C.L. 

It is not excluded that the data contain a mixture of 
two heavy quarks, one with \q\ = 2e/3 and one with \q\ = 
4e/3. We perform an unbinned maximum likelihood fit 
to the observed set of qi in data to determine the fraction 
p of exotic quark pairs. The likelihood of the observed 
set of qi can be expressed as a function of p by 

-^data 

L (p, q)= 1] 0- - P) p su(qd + pP cx (ft) (3) 
i=i 



Systematic 


Observed 


Expected 


Statistical uncertainty only 


95.8 


95.3 


+ Fraction of cc events 


95.8 


95.2 


+ Charge-flipping processes 


95.7 


95.2 


+ Weighting w.r.t. pr and y spectra 


94.4 


94.1 


+ Fraction of flavor creation 


93.7 


93.4 


+ Statistical error on Pf 


93.3 


93.1 


+ Jet energy calibration" 


92.4 


91.8 


+ Top quark mass 


92.2 


91.2 


"Reference |2'J. 



TABLE I: Expected and observed confidence levels as func- 
tion of the cumulated systematic uncertainties. 



Figure 1 (c) shows — In L as function of p. We fit p — 
-0.13 ± 0.66(stat) ± O.ll(syst), consistent with the SM. 
Using a Bayesian prior equal to one in the physically 
allowed region < p < 1 and zero otherwise, we obtain 
< p < 0.52 at the 68% C.L. and < p < 0.80 at the 
90% C.L. 

In summary, we present the first experimental discrim- 
ination between the 2e/3 and 4e/3 top quark electric 
charge scenarios. The observed top quark charge is con- 
sistent with the SM prediction. The hypothesis that only 
an exotic quark with charge \q\ ~ 4e/3 is produced has 
been excluded at the 92% C.L. We also place an upper 
limit of 0.80 at the 90% C.L. on the fraction of exotic 
quark pairs in the double tagged lepton-plus-jets sample. 
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